CLNS 96/1400 

RIMS-1063 

hep-th/9602074 



Multi-point Local Height Probabilities 
in the Integrable RSOS Model 

Sergei Lukyanov^'^'"^* and Yaroslav Pugai^''^ ** 

^ RIMS, Kyoto University, Sakyo-ku, Kyoto 606, Japan 
2 Cornell University, Ithaca, NY 14853-5001, USA 

and 

^ L.D. Landau Institute for Theoretical Physics 
Chernogolovka, 142432, Russia 



Abstract 

By using the bosonization technique, we derive an integral representation for multi- 
point Local Hight Probabilities for the Andrews- Baxter-Forrester model in the regime III. 
We argue that the dynamical symmetry of the model is provided by the deformed Virasoro 
algebra. 



February, 96 



* e-mail address: sergei@hepth.cornell.edu 
** e-mail address: slava@itp.ac.ru 



1. Introduction 



We start with an extended formulation of the problem which is attacked by solving 
in this work. 

1.1. General RS OS model 

Knowing the equilibrium structure of the crystal- vapor interface is a prerequisite for 
understanding and controlling the growth of crystals from vapor or dilute solution. At 
low temperatures the interface is expected to be atomistically fiat to minimize the surface 
energy, whereas at high temperatures it is rough and contains steps producing the config- 
uration entropy. This change in the structure is the roughening transition of the interface. 
The simple model generally accepted to depict the necessary minimum of configurational 
degrees of freedom describing the roughening transition is so-called "Restricted solid-on- 
solid" (RSOS) model 0. In this model crystal atoms are supposed to be located at the 
square lattice sites numerated by variables a without allowing an empty space in the bulk 
of the crystal. The interface is characterized by the height or number of crystal layers 
rria measured from some reference level. The statistics of the model is governed by local 
positively defined Boltzmann weights 
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assigned to every configuration {maTmbTmc^md) of heights round a face with the sites 
(a, b, c, d) ordered clockwise from the upper left. The restriction condition on the weights 
( p. . 1| ) expresses that no strongly fluctuating configurations are allowed. More precisely, if 
m's take integer values between 1 and r — 1 , 

1 < ma < r - 1 , (1.2) 

with the integer r > 4, then (|l7l|) is non zero only if any two neighboring heights around 
the face differ by ±1 : 

\ma - mb\ = 1 . (1.3) 
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The Boltzmann weight of a given configuration is a product of weights assigned to each face 
of the lattice. At low temperatures the interface tends to take on a height configuration 
minimizing the total surface energy 

£= e{ma,mb,mc,md) . (1.4) 

faces(a,b,c,d) 

Such configuration is usually referred as a ground state. In principle, the system may leave 
a room for a number of different ground states and we enumerate them by an integer 1. 
Consider the following lattice. 




The bullets here denote boundary heights which would be fixed somehow. In the present 
work we assume that the boundary heights take the same values as those in one of the 
ground states 0. So, there is one-to-one correspondence between these boundary condi- 
tions and the ground states and they can be labeled by the same index /. A partition 
function for the RSOS model with the designated l-th boundary condition is given by 
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Here is a number of the lattice faces. Of particular interest is the partition function 
Zi^^ at the thermodynamic limit N oo, when the lattice uniformly increases in all 
directions. The short distance interaction ensures, that 

^ K^Zi , A ^ oo , (1.6) 

and K is a partition function per face. The above mentioned boundary conditions break 
explicitly the translation invariance, hence Zi in ( p..6|) may depend on them, in contrary 



^ The translation invariance can be restored by a proper averaging over all possible values of 
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Other important quantities for studying are n-point Local Height Probabihties (LHP). 
They can be defined as foUows; Let us specify n reference sites belonging the same vertical 
column and label them from down to up by 1, 2, ...n. The probability that heights at these 
sites have the values 1 < /ci, < r — 1 reads: 



heights^my J = l 



n w 

faces(a^b,c,d) 



Again, of prime importance are LHP at the thermodynamic limit 
1.2. ABF model 
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The Boltzmann weights (|1 . 1|) are free parameters for the RSOS model, and an exact 
analytical calculation of LHP ( |1.8| ) for their general values does not seem to be a realizable 
problem nowadays. Nevertheless, there are some special cases when the model is a solvable 
one. As it usually is, a local integrability condition is given by the set of algebraic equations 
for the Boltzmann weights 0, [SI 
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Here u is some variable parameterizing the Boltzmann weights. Andrews et al. have 
succeeded in finding of solutions of the Yang-Baxter equation (YBE) (|1.9| ) for any integer 
r > 4. In addition to the spectral parameter u, the ABF solution depends on a parameter 
p. It gives a two parametric family in a manifold of the Boltzmann weights of the RSOS 
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model. The integrable weights read exphcitly: 
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Here we use the short notation 
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([mil][m-l])^ [w] 



(1.10) 
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[w] = /i(w)//i(l) 



(1.11) 



with 



+ 00 



h{u) = 2 \p\^/^ sm{nu/r) JJ (l -p") (l - 2p'^cos(27rw/r) + p^"") 



n=l 



In order to provide the convergence of the infinite product, the parameter p has to be 
restricted to the domain 

-1 < p < 1 . 

The positive factor R in ( |1.10|) determines a normahzation of the weights. 

The ABF weights family consists of two distinct manifolds with < m < 1, and 
— 1 < u < 0. According to 0], the manifolds are divided also into phases (regimes) by the 
line of critical points p — 0. In this work we restrict our attention to the so-called regime 
III: 

0<p<l, 0<w<l. (1.12) 

The parameter p measures a deviation from the criticality i. The value p = 1 can be 
associated with the zero temperature, since in this case the system is frozen in one of the 
ground states. As it follows from ( |1.1C1| ), each ground state in the regime III has the form: 
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^ In the case of the regime III, the scahng behavior in the vicinity of the critical hne is described 
by the Quantum Field Theory and referred as the Restricted Sine-Gordon model H. 
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Therefore, we identify the integer number /, enumerating the ground states, with min{li, I2) 
and 

l<l<r-2. (1.13) 
The same number specifies also the boundary conditions, used in the definition of LHP 

Now, let us break down the n — 1 rows between the reference points a = 1, 2, n 
and replace the spectral parameter u in these rows to < wi, w^-i < 1, respectively. 
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As well as in the homogeneous case, one can introduce n-point LHP Pki,...k„{ui, Un-i\l) 
at the thermodynamic limit on the lattice with such kind of dislocations. In this work we 
are going to study these quantities. Our main result is an integral representation for 
them given by the formulas (5.11). Notice, that it would be convenient for us to use the 
parameters ( and x 

^ 1 . 
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(1.14) 



rather then u and p. Hence, the n-point LHP for the inhomogeneous lattice are denoted 
as 

-Pfci,...A;„(Cl7 Cn-l|0 

in what follows. It is worth to point out here, that due to YBE ( |l.y| ) these functions 
depend on the ratios CjCi^ only 0. 



2. Corner Transfer Matrix and one-point LHP 

2.1. Corner Transfer Matrix 

The key calculation tool for LHP is the Corner Transfer Matrix (CTM) 0, [g, @]. 
In this section we concentrate on finding the one-point function to give an idea on this 
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important notion. Consider the square lattice with the reference site O as being the 
central one. 
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The probability Pkil) that the central height equals to k is determined by (pr7|), ( PT^ ) 
with n = 1. One can divide the whole lattice into four quadrants by the horizontal and 
vertical lines intersecting at O. The idea is to calculate the sum in ( |1.7|) in two steps. First, 
let us perform summation over the sites inside the quadrants with given values of heights 
along the cuts. Doing so, one treats the resulting expressions as matrix elements of the 
four Corner Transfer Matrices A, B, C, D associated with the corresponding quadrants. 
The second step is to find the sum of heights on the boundaries between the quadrants. 
More explicitly, consider the lower right quadrant with the hight configuration 

m = (mi, 7712, ...^1, /2) , mi = k 
along the horizontal boundary and 

for the vertical one. 

k m'^ m'^ m'^ . . . Ii 1 2 
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The Corner Transfer Matrix A(Q is a matrix with elements [A] ^ given by the partition 
functions of the quadrant with the fixed boundary heights. Notice, that we specify the 
dependence on the spectral parameter chosen in the form ( |1.14|) . CTMs B(C), C(C), D(C) 
for other quadrants are defined similarly. By the definition, CTM breaks up into r — 1 
diagonal blocks with the given value of the central height. Let P/c be the diagonal matrix 
which diagonal entries are units for the block with mi = k and all other elements being 



zero. Then the remaining sum in (|1.7|) over the heights at the cuts consists in taking the 
trace: 



Pk{l)=Tr 



Pt DCBA 
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(2.1) 



This expression does not depend on a normalization of CTMs. 

According to CTMs and their products would be well defined operators at the 

thermodynamic limit, if the Boltzmann weights are normalized in such a way that the 
partition function per face k ( p..6|) equals to one. For the ABF weights ( p..lO|) in the regime 
III this normalization means [§, [Q: 
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where the notation 
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{fci}=0 



is used. The function -R(C) obeys the set of relations 

R{C)R{C^) = 1 , 
R{x'C') = {xC-T 



With this normalization, the ABF weights satisfy 
- unitarity relation 
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(2.5) 
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"initial" condition 



W 



ki k2 
k4 ks 



C = l 



(2.6) 



From this point on we assume the normahzation ( p.2|) holds and at the thermody- 
namic hmit CTMs are weU defined hnear operators acting in the space covered by vectors 
(/c, 7712, ma, .../i, /2, h, ...)• We denote such spaces as Ci^k- 
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A(C) : 

2.2. Spectrum of CTM in ABF model 

At the thermodynamic hmit the Corner Transfer Matrix for solvable statistical models 
has the remarkably simple form 0,12] 



A(C) = C"^ 



[2.1 



where the Corner Hamiltonian He is a ^-independent operator. Notice, that this is a 
common property, since it follows from YBE (|1.9| ) and the crossing symmetry condition 
for the Boltzmann weights. The relation (|2.5| ) allows one to express the CTMs B, C, D 
in terms of A: 



B(C) = S,,fc A(x2c-i) , C(C)=A(C), D(C) = S;,, A(a;\ 



2/— 1^ 



(2.9) 



The matrix Si^k depends on the value of the central height k and the boundary condition 
only: 

^Lk = 



[k][l]I 



(2.10) 



Here I is the identity matrix. An explicit form of the function [/] is inessential for the 
problem under consideration, since it is cancelled out in the final expressions. 
Using the formulas (|2.8|) - (|2.1C1| ), we can rewrite ( |2.1|) as: 

r-l 



P,(/) = [k] Trc,^, 
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(2.11) 



Therefore, the calculation is reduced to finding the spectrum of the Corner Hamiltonian 
He in the space Ci^k- At the thermodynamic limit the spectrum is notable for the following 
reasons; It is 

— hounded from below , 

— discrete , (2-12) 
— equidistant . 
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These features seem to be rather general ones for solvable lattice models 0, [Qi. In the 
ABF case, they follow from the fact that the Boltzmann weights are quasiperiodic functions 



under the changing C e ( . Andrews et al. |^ have shown that the eigenvalues of the 
Corner Hamiltonian He in the regime III have the form: 



Az,fc - — + m 
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(2.13) 



(2.14) 



4r(r — 1) (r — l)r 

The discreteness of the spectrum implies that it does not change discontinuously with 
the temperature parameter p. Hence, to find the multiplicities of eigenvalues ( |2.13| ) 
of the Corner Hamiltonian acting in the space Ci^k ( ^-131 ) it is enough to perform the 
computation for any value of p ( |1.12| ). They can be derived at low temperatures limit 
p —>■ 1 hj combinatorial methods 0]. The amazing observation [0, fill], fl^ is the 
spectrum of He coincides with one of the Virasoro algebra generator Lq — ^ acting in 
the irreducible representation with the highest weight Ai^k and the central charge c ( |2.14|) . 
This immediately leads to the formula for the one-point LHP B 
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Here we denote the character of the irreducible module [jT^ by 
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and use the notation for the elliptic function 



Eiz, q) = (<?; q)ooiz; q)ooiqz ^;q)c 



(2.18) 



^ It is a starting point for so-called Angular Quantization approach for the integrable Quantum 
Field Theory [Q. In the scaling hmit the Corner Hamiltonian coincides (up to a multiplicative 
normalization) with the Lorentz boost generator, acting in the Angular Quantization space. 

^ Notice, that Zi = [I] Zi, where the function Zi is defined by (1.6) and [I] is the same as in 



ill). 
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3. Vertex operators and multi-point Local Height Probabilities 

3.1. Vertex operators 

Now we begin to study the multi-point LHP. In this subsection we follow the ideas of 
the work Consider the inhomogeneous lattice. Our challenge is to find the probabili- 
ties that heights in n successive sites from the same vertical column take values ki, . . .kn- 
It is useful to divide the lattice into 2n + 2 parts, as it is shown in the picture: 
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Again, at the first step of the calculation we specify the heights along the cuts and introduce 
the partition functions for the corresponding parts. The four of them are CTMs. Others 
are given by the products of the weights in the form of half infinite lanes: 



[ ^ i C )] 
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We treat such half infinite products with a given values of the spectral parameter (, as 
matrix elements of operators ^^(C) 5 The operators intertwine the spaces Ci^k 
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(3.1) 



and £z,fc±i 

and are referred as vertex operators ( VO) . The second step is to carry out the sum of heights 
over the cuts assuming that variables at the sites 1, ...n equal to /ci, ...kn , respectively. 
As before, it is equivalent to taking a trace. In such a way, we obtain the expression for 
the multi-point LHP: 



r-l 



i^.„...fc„(Ci,-Cn-ilO= E ^^^^"^ D(C)C(C)B(C)A(C) 
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(3.2) 
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D(c)**'^HCi)-**""-nCn-i)C(C)B(C)*'^"-nCn-i)..*"^(Ci)A(C) 



where as = kg+i — kg. 

To proceed further, one need to study the properties of VO. They can be derived from 
the heuristic graphic arguments . 



/. Commutation relations 

This property follows from the Yang-Baxter equation ( |1.9|) . 
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The face added in the first figure is pushed up to the infinity due to YBE. At the operator 
notations this implies that ^^(C) constitute the Zamolodchikov-Faddeev algebra in the 
Interaction Round a Face form: 



^kA—ko 



(Cl)*'=-'"(C2)|£«, 
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(3.3) 



//. Homogeneity condition 

Consider the composition of CTM and VO: 
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At the thermodynamic hmit this would be CTM A.{Q defined, however, in the space 
Ci^k±i rather than in ^ . To obtain an equahty, these spaces must be intertwined. It 
is provided by action of the operator $^(1) on the space because of the "initial" 

condition ( |2.6| ). Therefore, we have 

*^(C) A(C)|r,,. =A(0 . 



This leads to the equation: 



(3.4) 



///. Normalization condition 

To obtain the last property of VO we need to use the unitarity ( |2.4| ) and crossing symmetry 



(2.5) relations. They allow one to get so-called normalization condition: 



^ =V^I . (3.5) 



Now, let us express CTMs B, C, D and VO via A and respectively, by 

the formulas ( |2.9|) and 

**±(C) = *^(C) . (3.6) 



The latter follows from the crossing symmetry condition (|2.5|) . Taking (|3.4|) into account 



it is easy to find that LHPs can be written in the form |14]: 
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(3.7) 



where as = ks+i — kg, and Zi is given by ( p.l6|) . 
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3.2. Two-point LHP 

The simplest example of LHP (|3.7|) with n = 1 was treated in the previous section. 
Now we would like to show that from the properties of VO and CTM we can deduce the 
two-point function. As well as the one-point function, the two-point LHP do not depend 
on the spectral parameter. Then, according to (|3.7| ), the probabilities that heights at two 
adjacent sites 1 and 2 are equal to /c. A; ± 1, respectively, are given by: 



X 



Pk,k±i{l) = Zf' V[k][k±l] Trc 
It follows immediately from ( |3.5D that these functions obey the equation i 

In addition, the two-point LHP has to satisfy the evident requirements: 

Pk,k-i{l) = Pk-i,k{l) 1 
Pi,o{l) = . 

The formulae (|3.9| ) and (|3.10|) specify the two-point functions uniquely: 

k 

Pk+iAO = Pk,k+i{l) = ^(-1)^— [m] . 

771 = 1 

This equation will lend support to the validity of our general result for the multi-point 
LHP . 



(3.8) 

(3.9) 

(3.10) 
(3.11) 

(3.12) 



4. Bosonization of Vertex Operator algebra 

In this way, the multi-point LHP are expressed in terms of the traces of VO. From the 
mathematical point of view, the vertex operators and the Corner Hamiltonian constitute 
the algebra (|3.3|) , (|3.4|) , (|3.5| ). It acts in the direct sum of the spaces Ci^k- Unfortunately, 
our "definitions" were rather heuristic than rigorous ones. Due to Andrews et al. 0], 
we surely know that Ci^k has a structure of a Z-graded space provided by the Corner 
Hamiltonian. The features ( |2.12| ) are characteristic ones for a spectrum of the grade 
operator in a Fock space. Our main idea is to construct a representation of ( |3.3D , (|3.4|), 
( p?5| ) in a direct sum of Fock spaces. This representation is expected to be reducible. 
However, its restriction to an irreducible component would provide the proper multiplicities 
for the eigenvalues of He. A similar bosonization procedure was applied for calculations 
of conformal blocks in the Conformal Field Theory [|l^ and lattice correlation functions 
for the XXZ Heisenberg spin chain |jT^ , |jl^ . 



Note, that (3.9) obviously follows from the definition of the one-point function and ( ^.15| ). 
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4.1. Bosonic representation of VO 

It is convenient for us to perform the simple transformation of VO 

*^(0I/:.. =^'-^ V^*^(C)k,. • (4.1) 

The operators also constitute the associative quadratic algebra. It can be easily 

verified by substitution of ( [4.1| ) into (|3.3|). In particular, the commutation relation of 
*+(Ci) and *+(C2) has the form: 

*+(Cl)*+(C2) = i?(ClC2"') *+(C2)*+(Cl) , (4.2) 

where the function R{C) is given by (|2.2D . Our first task is to carry out the bosonization of 
^"•"(C) satisfying this simple equation. To do it we follow the procedure developed in 0; 
Let us introduce the operator valued function ip{C) obeying the commutation relation: 

[y'(Ci),V'(C2)] = -lni?(CiC2"') • (4.3) 

The quasi-periodicity of -R(C) under the substitution C e^'^'^Q implies the following 
Laurent decomposition for ^{Q) '■ 



^{0 = -f—l{Q-^V K) +^Y,—C'-■ (4.4) 



In order to provide (|4.3|) , we specify the non- vanishing commutation relations for the 
operators "P, Q, as being 

H^[(r-l)m]^ 

[2m] ^[rm]^ (4.5) 

where the notation 



is used i. With these definitions the operator v'(C) satisfies the commutation relations 
It is easily checked, since the function -R(C) (|2.2|) can be rewritten in the form: 



It should not be confused with the symbol 
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The Heisenberg algebra (O) is represented in Fock spaces. In the usual fashion, one 
defines the Fock space Tp as the module generated by action of the creation operators 
fi-ni n > on the highest weight vector vp 

/3„vp = 0,n>0; Pvp = Pvp. (4.7) 

Namely, Tp is covered by vectors /3_„^/3_„2 ■ ■ -(^-nj^p , > . It is endowed with 
the structure of a Z-graded module by the operator: 

E[2ml T [rml T 1 , , 

Nx[(r-l)m]^ 2 24 

Its eigenvalues in jFp have the form ( ^.13| ) and we would like to identify ( |4.8D with the 
Corner Hamiltonian. As it follows from ( |2.13p , only the Fock spaces 

^l,k = J^ Tl-{r-l)k , (4.9) 
V'2(''-l)'' 

with / and k being integers, are relevant in our construction. 
Now, let us introduce the following bosonic representation: 

*+(C)=e^^(^): J-^,, (4.10) 

Then the relation ( [4.2| ) is satisfied immediately. In order to describe the bosonic repre- 
sentation for the vertex operator ^'"(C), we need the notations: 

[/(C) = e*'^(^) , U{z) = e-^M^^)+'^(^^"')) . (4.II) 
Define also the operator F{z,V) depending on the "zero" mode V 

F(z, V) = ,— ^^r.i.+i)-. ) ^412) 



where uj = y ^^7-^ and E{z, q) is given by ( p.l8|) . We propose the following bosoniza- 
tion prescription: 

*-(C) =r'i4^ U{OU{z) F{zC\V) = / F{Cz-\ -V) U{z)U{0 ■ 

Jq ZTXIZ Jq ITIIZ 

(4.13) 

Here the anti-clockwise contour C encloses the poles z = C,x^^'^'^'^ (m = 0, 1,2, ..) of the 
integrand and 77 is some constant, which we are going to specify later. It is significant 
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that the integrand in ( ^4.13| ) is a single-valued function of z on the complex plane, so the 
integration contour is closed. Therefore, the action of the operator (^4.13|) is well defined 
for an arbitrary Fock space Tp. We claim that; 
The operators ^4AW , 



(4.14) 



satisfy the commutation relations 

*°(Cl)*'(C2)^„ = 



c-\-d=a-\-b 



k + a + b k + c 
k + b k 



C1C2"' 



*'(C2)*^(Ci)l^,,. • (4.15) 



The coefficients U ( |4.15| ) are connected with the Boltzmann weights (|1.10| ) via the simple 
transformation 
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m4 ma 
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w 



mi m2 
m4 ma 



(4.16) 



It is apparent that the commutation relations ( f4.15| ) are equivalent to (^]5) if VO 
and are related as in (|4.1| ). The explicit form of the matrix ( [4.16| ) is presented in 
Appendix A. We give the flavour of the proof (|4.15|) in Appendix B. It is easy to see that 
the operators (|4.10| ), (|4.13| ) obey the proper commutation relation ( |3.4| ) with the Corner 
Hamiltonian (|4.8|). One can also check that they satisfy the normalization condition: 



I *-(a:-2C)*+(C)-*+(x-2C)*-(C)}|^^^ = A (-1)^-^ [k] I 



(4.17) 



where the numerical constant A is equal to 



2r-2 ^2r^ 



(4.18) 



The occurrence of the multiple ^Jl — z\z^i is connected with a fact that we do not use the 
normal ordered forms of the exponents (|4.11|) . The normal ordering of U{Q produces 
only an additional multiple factor which is a non zero and finite constant for < x < 1. 



At the same time, the normal ordering of U{z) gives rise to the factor a/1 — z\z^i- It can 
be cancelled out by a suitable choice of 77. We specify this constant in such a way that 



A 



and ( ^l.lTp would be equivalent to the normalization condition (|3.5|). 
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4-2. Felder complex 

Now we should concentrate our attention on the description of Ci^k in terms of the 
bosonic Fock spaces. RecaU, the Corner Hamiltonian was identified with the grade operator 
( [4.8| ). Although the eigenvalues of ( [4 .81 ) coincide with ones for the Corner Hamiltonian, 
their multiplicities are different. Hence, the spaces JF; ^ and Ci^k cannot be identified 
directly. More sophisticated treatment shows that the bosonic representation of the VO 
algebra is reducible. To construct the proper space of states one has to extract an ir- 
reducible component by throwing out some states from the Fock spaces. Explicitly, the 
procedure reads as follows; Introduce the notation 

V[z) = 6"*'^'^*^^^)+'^*^^^"^)) (4.19) 

with given by 



^^0-J^^{Q-'Vln()-,j:^C- (4.20) 



m/0 



and 



[(r - l)m]r, ara = [rm]^ [5^ . (4.21) 
A central place in the analysis belongs to the operator 

X = r' (f^Viz)F{z,V), (4.22) 
J Zmz 

where 

rr ( ^-l-2ir-l)uj ^ ^2r-2\ 



with u; = — Y ;~rr^ and the constant ry provides the regularization of V{z). The operator 
( [1.22|) defines the following maps: 



X2j — '■ ^l-2j{r-l),k ^-l-2j{r-l),k , 



(4.23) 

-yr — l — l T T" 

2j+l — ■ J--l-2j(r-l),k ■^l-2{j + l){r-l),k 



for jGZ, l<^<r — 2, 1 < k < r — 1. As a result, we can construct the infinite chain: 

^ ^ . (4.24) 

^ J'2r-2-Lk ^ -ri.k ^ -r-l.k ^ ••• 



We claim that it is the Felder resolution fl8|. In other words 
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/. The chain of maps U-^A) a BRST complex for any < x < 1 , i.e. 



XjXj_i = ; (4.25) 
//. The cohomologies of the complex turn out to he non-trivial only for j = 

Ker Xj I Im Xj_i = , if j ^ . (4.26) 



The first statement can be checked by the direct calculations []T^. The arguments in 
favour of ( [4.26|) are based on the fact that the operator ([4.22| ) commute with He ([4.8|). 
So, the cohomology spaces are Z-graded, as well as jFp. Let us consider the finite- 
dimensional eigensubspaces of the operator He in Ker Xj I Im Xj-x. The restriction of 
( [4.24|) on these subspaces provides BRST complexes involving only a finite number of the 
maps. Dimensions of the cohomology spaces for such complexes are integers and would 
not depend on the continuous parameter x. Therefore, it is enough to find them for an 
arbitrary < x < 1. The statement ( [4.26| ) surely holds at the limit x — 1, where our 



construction turns out to be the Feigin-Fuks-Dotsenko-Fateev bosonization [^^, |jT^. In 
addition, we conclude that the spectrum of H^ in the space 

= Ker Xo / Im X_i = Ker^, , X' j /m^,_,_,,, X^-^-' (4.27) 

coincides with the spectrum of the grade operator Lq — ^ in the corresponding irreducible 
representation of the Virasoro algebra. 

To complete the identification of the factor-space ( [4.27|) with one notes that VO 
satisfy the commutation relations: 

*±(C)X,=X, *±(C). (4.28) 

So they are BRST invariant operators. 



The BRST properties and the structure of the complex ( 4.24 ) lead to the following 
computation procedure of traces over Li^^ in terms of traces over the Fock spaces |T^. 
Introduce the notation: 

Or. = *-'^n^"'Cl)-*"""-n^"'Cn-l) *""-nCn-l)...*"nCl) • (4-29) 
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Then, due to the commutation relations ( f4.28| ) the foUowing diagram is a commutative 
one: 



X_2 X_i 
^ ^-l+2r-2,k ^ 

> ^-l+2r-2,k ^ 

X-2 X-i 



Lk 



i,k 



^-i,k 

^-l,k 



Xn 



(4.30) 

The complex is exact excluding the j = term, hence the trace of the operator x^^'^On 
over the space ( [4.27| ) equals to the following alternating sum: 



Tr 



2r-2-i,fe 



!-2j(r-l),fe 



!~2j{r-l),fe 



(4.31) 



j = -oo 



j = -oo 



In the simplest case 0\ = \ the formula (|4.31|) gives 



„4Hr 



!-2j(r-l),fe 



j — — oo j — — oo 

(ri-(r-l)fc+2(r-l)rj)2 

j = — oo j= — oo 

It coincides with the character xz,fc(x^) ( |2.17| ), due to the Jacobi identity: 



„4Hf 



X 



(rl + (j l)k + 2(r-l)rj)^ 

/jQ (r — 1 ) r 



(4.32) 



+ 00 



(m — 1 ) m 



m= — oo 



5. Integral representation for multi-point LHP 

5.1. Calculation of traces 

In this section we work out an integral representation for LHP by using the results of 
the bosonization procedure. In terms of VO ^^(^) the equation ( p.7|) can be rewritten 

as 

n-l 



s=2 



(5.1) 



Trr,,^ x^"^ *-'^i(x-^Ci)-*"""-^(a;-'Cn-i)*""-^(Cn-i)...*"^(Ci: 



-o-i 1-^-2 
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It is obvious that in the calculation of (|5.1|) we can perform the traces over the zero and 
oscillator modes separately. The zero mode contribution is governed by the structure of 
the complex ( [4.3U| ) and it was obtained in |T^ . The trace over the oscillator part is carried 



out by using the Clavelli- Shapiro technique [jT^. More explicitly, the prescription 



implies the introducing together with /3„ the oscillators 

^*m-- [/3^,/3:] = o, 

satisfying the same commutation relations as ( |4.5| ). Define the following operators acting 
in the tensor product T[l3\ ® of the Fock spaces: 

h^ = {l- x^"^)-^ ® 1 + 1 ® /?* ^ , m > ; 

(5.2) 

6^ = ® 1 + (x^"^ - 1 ® /3* ^ , m < . 

Notice, that we omit the index P since the result does not depend on the zero mode 
eigenvalue in this case. For any operator 0[l3\ acting in one can construct the 

operator 

0[h] : T[[3]®T[[3*] ^ T[f3]®r[P*] 

by the replacement of the modes fim with hm (|5.2|) . Then the "oscillator" trace over 
the Fock space is expressed in terms of the vacuum expectation value << 0\b] »=< 
0\O[b]\0 > of the operator 0[b] with respect to the vacuum |0 >= v v. Namely, 



Tr 



osc 



{x\x')-J «0[h]» . (5.3) 



Due to the Wick theorem, the expectation value of a product of the exponential operators 
( [4.11|) is factorized into the two point functions @: 



«[/(C2)f/(Ci) »=7^2 G'(CiC2"') 

« u{C2)u{Ci) »=vnn ^^'(CiCs"') (5-4) 

« f/(C2)f/(Ci) »= ^'(CiCs"') • 



The calculations of the functions and constants in (O) are straightforward, and their 
explicit forms are written down in Appendix C. 



These functions were calculated in the work (see also [p2|| ), and we preserved here the 
same notations. The prime should not be confused with the symbol of the derivative. 



20 



5.2. Integral representation 

We begin with the following trace: 



(5.5) 



According to our analysis, ( |5.5| ) admits the bosonic representation ( |4.31D with 

dz 



27riz 



F{C2Z-\-V) U{z)U{C2)U{(:i) . 



(5.6) 



First, let us derive the zero modes contribution into the integrand. For the j = term in 
the sum ( [4.31|) one obtains: 



2i^/^{Q-iV\nz) -i^/^(Q-iPlnC2) f^-i^/^iQ-iV In Ci) 



^ r— 1 6 



(5.7) 

where u> = P = I — ^—^k. All other terms can be found by the substitution 

/ ±Z — 2j(r — 1) , J G Z into ( |5.7|) . Then, performing the sum over j we get the zero 
mode contribution: 



where 



[xSx^U (C1C2"') ' 0;,fc(ClC2^-') «2"')""^ F{C2Z-\-P) 



A,fc+i(CiC2"';0 = ^r' ^ (-1)^-'-' ^ G'(CiC2"') (C1C2"') 



1 \ 4r 2 



X 



^l = IC2 



2'niz 



(5.8) 



The trace over the oscillator modes is provided by the Wick theorem and ( ^.4|) . As a 
result, we derive the integral representation for ( p.5| ): 



0,,fc(ClC2^-') I^'(Cl^-')Vr'(C2^-^) «2"')'^-'^ F{C2Z-\-P) . ^^'^^ 



The explicit form of the functions and constants appearing here is rather complicated 



(see Appendix C). The function (|5.9|) represents LHP on the lattice with more general 
dislocation, than one considered before: 
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k 












• 


k 








B 





















In the case C,\ = x^C2 the expression ( |5.9| ) is simphfied drastically. Indeed, using the 
obvious relation 

it is possible to show that the two-point LHP (|3.8|) are given by: 



( 1 1\6 



X 



r _dz_ ^ 2 E(^-lx^'^'-^("-^)^a;'") 



(5.10) 



In Appendix D, we prove that these functions coincide with ( |3.12| ) 



In a similar manner, we have obtained the integral representation for the multi-point 
Local Height Probabilities: 

kji — k-i -\-n — 1 — 

Pfc„...fc„(Ci,-Cn-iio = (-1) — ^ — [ki] zr^x 

n—l n— 1 „ , n— 1 

n «i n n { £ ^} ( n c^.-) x 

n-l 

s<j i<s j>s s=l 



where 



^ ^ E{z,x^) ' 
' E(z,a;2) 

^(^-2H+2(r-l)fc^^2r) 
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and the functions [/c], Zi, E{z,q), Oi^kiC) are given by (|LTT|), ( pJ6D , (pJSD , 
respectively. The integration contour Cg, {s = 1, ...n — 1) in ( |5.11[ ) goes anti-clockwise and 
encloses the following poles: 

^,^r{2m+l) + {r-2}{k,-kj + i) ^ ^ ^_^.^r (2m+l) + (r-2) (A;,+i - fc, ) y ^ 

where m = 0, 1, 2, ... . 



6. Conclusion 

In summary, we would like to make some comments on an algebraic interpretation 
of the presented results. Our construction leads to a natural conjecture that the space 
Ci^k can be treated as the irreducible representation of some infinite-dimensional algebra. 
Since at the limit x — 1 the bosonization turns out to be the Feigin-Fuks-Dotsenko-Fateev 
procedure then we have called the algebra as the deformed Virasoro one ViVc^x 

in the work [^. The explicit basis for Vir^^x was found in the works [^5|i. It can 
be described in terms of the oscillators A^, : 



[(r - l)'m]x otra = [rm]x Pm = (x - x ) m X 



(6.1) 



If one introduces the field 



A(^) = xV2-(-i)7' 



exp I - ^ X^z~ 



(6.2) 



then the generating function for the elements of Vir^x has the form ||2^, [|2J], [|25| 



T{z) = A{zx-^) + A-\zx) . 

This field commutes with the operators ( [1.23| ) and obeys the relation 

fiCz-') T{z)TiO - fizC') T{C)T{z) = 

2iT {x - X-') [r - l]x [r]x ( S{Cz-'x-^) - d{Cz-'x^) ) 

where we denote ^(-2) = 2^ Sm=-oo 



(6.3) 



(6.4) 



f{z) = il-zy 



„2-2r. 



( yrp2r+2. rp4:\ /'y^4— 2r.^4^ 



J. Shiraishi et al. |24] use the notations p = x , q = x ^ , t = qp 
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The limiting x ^ 1 behavior of ( |6.3| ) is governed by the expansion: 

T{z) = 2 + r(r - 1) (x - x-'f ( ^ L^z""^ _ ^ j + 0((x - x-^) , (6.5) 

\m= — oo / 

and (|6.4| ) gives us the Virasoro algebra commutation relations for Lm with the central 
charge c (|2.6|) . 

Remind, that the operators ( [4.10| ) , (|4.13| ) intertwine the irreducible represen- 



tations Ci^k with the same first index /. In the terminology of the work |jT^, they are 
the first kind VO. It would be rather natural to introduce the second kind VO, which 
intertwine the spaces Ci^k and £z±i,fc • Their bosonization is easily performed in terms 
of the oscillators a^i ( [4.21D . Then the " Corner" space of states i 



nz= Yl ®'^^'^ (6.6) 

l<l<k<r-l 

of the ABF model at the regime III can be treated as the irreducible representation of the 
complete algebra of both kinds VO. Notice, that we restrict the sum ( |6.6| ) by 1 < / < 
/c < r — 1, since 

^l,k — ^r — l — l,r — k • 

The integrability of the ABF model implies an existence of an infinite dimensional 
Abelian symmetry. The deformed Virasoro algebra provides another type of the symmetry 
which is usually referred as the dynamical one [1^] (see also ||2^). The generators of 



such symmetry do not mutually commute, but create the space of states of the model. 
Other well known examples of the dynamical symmetries were produced by the Conformal 
Field Theory [^. From the algebraic point of view the symmetry corresponding to the 
deformed Virasoro algebra is very similar to the conformal invariance, although the simple 
geometrical meaning of the latter seems to be lost. 
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The lattice analogue of the Angular Quantization space 
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7. Appendix A 

For convenience, we collect in this Appendix two different forms of the matrix U 
( [4.16| ). In terms of the parameters u and p, it reads: 



U 
U 

u 



m± 2 


mil 


mil 


m 


m 


mil 


mil 


m 


m 


mil 


m =1= 1 


m 



R 



R 



-R 



[m i u] 
[1 — u] [m] ' 
m] [m =1= 1] 



[1 — u] [m] 

Via the variables ( and x (|1.14 ), the same functions are given by the formulas: 



U 
U 
U 



mi2 


mil 


mil 


m 


m 


mil 


mil 


m 


mi2 


m =1= 1 


mil 


m 



( r — 1 ) ( =p TTi — 1 



) E{x 



2r-2 ^2r 



T2(r-l)m^^ ^2r 



2r\ 



^2(r-l)m ] 

-i? x^ — ? — C" 



^(^2(r-l)(±m+l)^^2r-) EjC^X^"-) 

The function = -R(C) is specified by the equation (|2] 



(7.1) 



(72) 



8. Appendix B 

We show that bosonic operators ( [4.1C1|) , ( [1.13| ) satisfy ( [1.15| ) by working out the ex- 
ample of the commutation relation 



*-(Ci)*+(C2)k,, =u 
u 



k k-1 
k + 1 k 



C1C2 



*+(C2)*-(Ci)Ia..+ 



k k + 1 
k + 1 k 



C1C2"' 



*-(C2)*+(Ci)|/: 



To simphfy the formulas, let us return to the parameterization (|1.14|) . In terms of the 
variables p and u the function F(C, V) = F[u, iv] ( [4.12| ) reads: 

[1/2 + u-ru;] 



2) 



where to 



[1/2 - u] ' 

^'^"^ V and the symbol [u] is defined by (|1 . 1 1| ) . Then, the following equation 



j.r , 11 [r-l][ruj + ui-U2][l/2 + v-Ui][l/2 + U2-v] 

F\v - Ui.ru - r + 1\ = — -+-^, -+-^, -x 

^ ^ [r-l + ui-U2][ruj][l/2-v + Ui][l/2-U2 + v] 

,11 [r-l-ru][ui-U2][l/2 + U2-v] 

F\v — U2,ru! — r + 1\ -— -— F\v — ui,u\ 

^ ^ [iv][r-l + ui-U2][l/2-U2 + v] ^ J 
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1.3) 



is a simple consequence of the famous Riemann identity: 

[2x] [2y] [2z\ [2t] = [x + y + z + t\[x - y - z + t\[x + y - z -t\[x - y + z -t\ + 
[—X + y + z + t\[x — y + z + t][x + y — z + t\[x + y + z — t] . 

By using ( |8.3|) , it is easy to obtain the relation between the exponential operators ( [4.11| ): 



(8.4) 



U{C^)U{z)F{zC^\V)U{C2)W,='^ 
F{zQ\V)\c,,,+\J 



k k + 1 
k + 1 k 



C1C2"' 



U{C2)U{Ci)U{z)x 



k k-1 
k + 1 k 



C1C2"' 



(8.5) 



U{C2)U{z)F{zC^\V)U{C^)\c, 



here ^1^2 

Now, to derive ( ^.1|) , we should integrate both parts of this 

equation over the variable z. 



9. Appendix C. Reference formulae for trace calculations 

Here we collect the formulae for the functions and constants in 



G'{0 



X 



(rf,2r/-. rf,2r ^4 ('t-4/"- 'r2r ^4 ^4^ 

('™2r+6/— 1. ™2r ™4 ™4^ /'™6/--l.™2r ™4 ™4^ 
('^2r+4/--l. ^2r ^4 ^4^ ('^8/--1.^2r ^4 ^4^ 

(rp2r-l/- 2r 4\ / 2r+3 /■- 1 . ™2r 4n 

S?*^ ^-^/cx)*^"^ S 7*^ 1 Jog 



( rtf . 2ir rir'4: \ ( ,t^5 /" 1 . /y. 2?^ ^j^4 ^ 

(C? )oO (*'' C 5*'' )oO 

(rr2r-2/-. ^2r^ (^l2H^27^^TTT3i2rV 

V'*' S 7 /OO S 7"^ A 



n 



I X 



7 loo 7 JOO 



V} = G"(l) . 
The functions satisfy the relations 

G'{Cx)G\Cx-^) = W'-\Q , 
W-'(Cx)W(Ca:-i) = G-i(C) 7 
G'(C) = G'{x\-^) , 

w^'(C) = vr'(x4c-') 7 

G'(0 = G'(x\-i) . 



(9.1) 

(9.2) 

(9.3) 

(9.4) 
(9.5) 



(9.6) 
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10. Appendix D. 

In this Appendix we demonstrate that the expressions ( |3.12| ) , ( ^.lOp for the two-point 
LHP are equivalent. The integral ( ^.101) is given by a sum of the residues of the poles 



situated inside the circle |2;| = x The direct calculation leads to the formula: 

i-k --1 

Pk,k+iil) = (-1) Zi — — J- — X 

E^x^"^ ^,x^'^)(x'^;x'*)oo (10.1) 

^rLo{LO-l) + ^ I j^^2H-2(r-l)A;j _ ^4rlk+2rl j ^^-2rl-2{r-l)k^^ 

Here A; jt, c are defined by ( |2.14| ) and 
J(m) = -^5^' (-1)^ x^^'-'^E{-m^x''^''-'>+^^'--'^'-\x^^'-'^')E'{mx-^%x^n . (10.2) 



s=0 



The prime in ( |10.2| ) denotes a derivative with respect to the first argument. Using simple 
properties of the elliptic function E{z, q), one can show that J(m) satisfies the identity: 

Jix^'^-^m) = -m-^£;(-m2x^('^-i)'^, a;^(^-i)^)E(m, x^"^) + m-^J{m) . (10.3) 

It provides the recursion relation 

Pk,k+i{l) + Pk,k-i{l) = {k\ XiAA , (10-4) 



where Xi,k{x'^) is given by (|2.17|) . It only remains to check that the function -Pi,o(0 obeys 



the condition ( |3.11| ). In terms of integrals this means: 

Eizx^^-^-^'^^.x^'^) 



I _ ^2r-2^4H+(4r-6)(r-l)^^8(r-l)r\ 

J^^^^^2 2Tciz ' E{z~^x^,x'^ 

f dz ^-2l^6ljp^ ^2r-2^-4H+(4r-6)(r-l) ^8(r-l)r\ E{zx'^'' ^''^ , X 

J\^\^^2 2Txiz ' E{z~^x'^,x'^ 



z\=x 



2r-3-2rl ^2r\ ^^^'^^ 



To prove the equality we should change the variable z ^ z ^ and then deform the 
integration contour at the left hand side of ( p.0.5|) . 
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